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Abstract 

The first generation of stars was formed from pri- 
mordial gas. Numerical simulations suggest that the 
first stars were predominantly very massive, with typ- 
ical masses M > 100M©. These stars were responsible 
for the reionization of the universe, the initial enrich- 
ment of the intergalactic medium with heavy elements, 
and other cosmological consequences. In this work, we 
study the structure of Zero Age Main Sequence stars for 
a wide mass and metallicity range and the evolution of 
100, 150, 200, 250 and 300M Q galactic and pregalac- 
tic Pop III very massive stars without mass loss, with 
metallicity Z = 10~ 6 and 10~ 9 , respectively. Using 
a stellar evolution code, a system of 10 equations to- 
gether with boundary conditions are solved simultane- 
ously. For the change of chemical composition, which 
determines the evolution of a star, a diffusion treat- 
ment for convection and semiconvection is used. A set 
of 30 nuclear reactions are solved simultaneously with 
the stellar structure and evolution equations. Several 
results on the main sequence, and during the hydro- 
gen and helium burning phases, are described. Low 
metallicity massive stars are hotter and more compact 
and luminous than their metal enriched counterparts. 
Due to their high temperatures, pregalactic stars ac- 
tivate sooner the triple alpha reaction self-producing 
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their own heavy elements. Both galactic and pregalac- 
tic stars are radiation pressure dominated and evolve 
below the Eddington luminosity limit with short life- 
times. The physical characteristics of the first stars 
have an important influence in predictions of the ioniz- 
ing photon yields from the first luminous objects; also 
they develop large convective cores with important he- 
lium core masses which are important for explosion cal- 
culations. 

Keywords first stars, stars: models, evolution 



1 Introduction 

A first generation of stars composed of primordial 
nearly pure H/He gas must have existed, since heavy 
elements can only be synthesized in the interior of the 
stars. These first stars, also called Population III (or 
Pop III), were responsible for the initial enrichment of 
the intergala c tic medium ( I GM) with heavy elements 
( Bond I fl98ll: iKlapp I Il98ll Il983l Il98i ICavrel I Il986 



19961: ICarr ||l987l 11994 iBromm et al. 112002}) 



UV photons radiated by the first stars, perhaps to- 
gether with an early population of quasars, are ex- 



(Haiman and Loeb Ill997: Ferrara 


19981 Miralda-Escude et al. hi 


2000l iTumlinson and Shull 1 2000: 


Brommetal. 2001: 


Schaerer l2002h. The enerev output from the first stars 



might have left a measurable imprint on the cosmic 
microwave background (CMB) on very small scales 



( Visniac 1 11987: Do delson and Jubas Ml995h . 



Howcver, despite an intense observational effort, the 
discovery of a true Pop III remains elusive. To probe 
the time when star formation first started entails ob- 
serving at very high redshifts z > 10. 

The first stars were typically m any times more mas- 
sive and luminous than the Sun (Larson and Bromm 
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20041 ) . A review of theoretical results on the formation 



of the first stars has been made by IjBromm and Larson 
20041 ). The masses of the first star-forming clumps 



would have been about 500 to 1,OOOM . Several nu- 
merical simulations suggest that the first stars were 
predominantly Very Massi ve Stars (VMS), with typ- 
ical masses M > lOOM^ dBromm et al. 1 1l999i 120021 



Nakamura and Umemura bOQlUAbel et al. I2000ll2002l )| 

and these stars had important effects on subsequent 
galaxy formation. 

In another scenario, the hypothesis that first stars 
were VMS (M > 140M Q ) has been strongly criticized 
because the pair-instability supernovae yield patterns 
are incompatible with the Fe-peak and r-process abun- 
dances found in very metal poor stars. Models includ- 
ing Type II supernova and/or hypernova from zero- 
metallicity progenitors with M = 8 — 40Mf^ c an better 



cal masses and Initial Mass Function (IMF)? and, 2. 
During their cuasi-static evolutionary phases, do they 
have radiation driven winds or mass loss due to other 
mechanisms? 

Star formation and accretion ca lculations suggest 
that t he first stars were very massive ( Omukai and Palla ^ 
120031) . On the other hand, by comparing the observed 
abundance patterns of Extremely Metal Poor (EMP) 
stars with supernova explosion calculations, some au- 
thors have concluded that the first stars are more likely 
to have masses in the ran ge ^ 20 — 30M^, but not 



explain the observed trends (jTumlinson et al. 1 12004). 



The same authors also pointed out that the sole gen- 
eration of VMS (M > 14OM ) cannot be possible and 
suggested that some VMS could be formed as compan- 
ions of stars with masses M < 140M Q . Their Initial 
Mass Function (IMF) proposition match quite well with 
the reionization and nucleosynthesis evidence. Tegmark 
et al. (1997) argued that the minimum baryonic mass 
is redshift dependent and lies in the range ~ 10 3 7 to 
10 6 M Q , for z ~ 10 and ~ 15, respectively, and that 
a participation of ~ 10~ 3 of the whole baryonic mat- 
ter in the generation of luminous stars is sufficient to 
reheating the universe by z ~ 30. 

The Wilkinson Microwave Anisotropy Probe has ob- 
served the large-angle polarization anisotropy of the 
CMB JCen 1120031; iKogut et al l liooj ISokasian et al. 
2003L IWvithe and Loeb 1 120031) . Some measurements 



have been interpreted as a signature of a substantial 
early activity of massive star (MS) formation at high 
rcdshifts z > 15. 

The supernova explosions that ended the lives of 
the first stars were responsible for the initial enrich- 
ment of the intergalactic medium with heavy elements 



l Ostriker and Gnedin Ill996 : Gnedin and Ostriker |[l997l | 



Bromm et al. 



200i 



j lYoshida et al. 112004) . An inter 



csting possibility unique to zero-metallicity massive 
stars is the complete disruption of their progenitors 
in pair-instability supernovae explosio ns, which are 
predicted to leave no remnant behind (IBarkat et al. 



1967 tlOber et al. ll983tlBond et al Hl984lFrver et al. 



200lt iHeeer and Woosleyl 120021 Irleger et al 



2003 



Bromm and Larson 1 120041 ). The later works consider 



that this peculiar explosion mode could have played 
an important role in quickly seeding the intergalactic 
medium with the first metals. 

Related to the first stars, there are two very im- 
portant unsolved questions: 1. What are their typi- 



more massive that 130M^ ( Umeda and Nomoto 112002 



Heger and Wooslev 112002 ) 



We have suggested that a possible solution to this 
inconsistency is that first stars are born very massive 
but that during their cuasi-static evolutionary phases, 
lose mass and reach the pre-supernova stage with the 
masses required from supernova ca lculations to repro- 



duce the EMP abundance pattern (jKlapp et al. 112005 



Bahena 1 12006). It is then very important to estimate 



the amount of first stars mass loss during the cuasi- 
static evolutionary phases. 

Kudritzki (2000, 2002) calculated wind models of 
massive stars between 100 and 300M Q and metallicities 
in the range 0.0001 to 1.0 solar, in an effective tempera- 
ture range from 40, 000 to 60, 000 K, with the objective 
of predicting mass-loss rates at very low metallicities 
applicable to the first generation of massive stars. It 
was found that for very low metallicities, the line driven 
mechanism becomes very inefficient and wind solutions 
can only be obtained very close to the Eddington limit. 
He also pointed out that very massive stars are pul- 
sationally unstable, which might contribute to stellar 
mass loss, in particular at low metallicity when the con- 
tribution of the radiative driving to the wind decreases. 
However, the critical mass for the onset of the nuclear 
pulsational instability is uncertain. 

Other mechanisms could induce first stars mass loss, 
for example, the low metallicity rotation models of 
Meynet and Maeder (2002) show fast rotating cores 
that lose significant amounts of mass and thus angu- 
lar momentum. Rotation by enhancing the luminosity 
and lowering the effective gravity increases the mass 
loss rate. Then, it is possible that pulsation, rotation 
and Luminous Blue Variables (LBV) type phenomenae 
could induce significant amounts of mass loss during 
the first stars cuasi-static evolutionary phases. 

Motivated by the above arguments, in a series of 
papers we will study the structure, evolution and nu- 
cleosynthesis of the first stars with and without mass 
loss and rotation. In this Paper I we present the evolu- 
tionary results without mass loss and with no-rotation. 

This work is organized as follows: In Sect. [2] we de- 
scribe the initial conditions of the stellar models and 
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the way in which the main physical variables are com- 
puted. Then, in Sect. [3] we describe the main results, 
and in Sect. 2]we discuss our results and compare them 
with other authors. Finally, in Sect. [5] we outline our 
conclusions. 



2 Numerical modelling and input physics 

We present Zero Age Main Sequence (ZAMS) models 
for stars in the mass range 1 — 10, 000 Af© with com- 
positions (X,Z) = (0.765, 1CT 2 ) for Pop I, {X,Z) = 
(0.749, 10~ 3 ) for Pop II, (X,Z) = (0.765, lO" 6 ) for 
galactic Pop III, and (X,Z) = (0.765, 10" 9 ) and 
(1.0, 10~ 10 ) for pregalactic Pop III. 

In this paper, evolutionary models for Pop III stars 
have been calculated without mass loss and with no- 
rotation. The chemical composition of the models is 
(X,Z) = (0.765, 10" 6 ) and (0.765, 10" 9 ), for galactic 
and pregalactic stars, respectively. For the evolution 
the chosen stellar masses were 100, 150, 200, 250 and 
3OOM . 

The main difference between galactic and pregalac- 
tic stars is their initial metallicity. According to Castel- 
lani (2000) we consider a Pop III metallicity range from 
Z = 10~ 6 to 10 -10 . The transition from pregalac- 
tic to galactic stars occur at the critical metallicity 



Z ~ 10~ 3 - 5 Z© (jBromm and Larson 1 120041) . Pregalac- 



tic stars are characterized by self-producing their own 
heavy elements; galactic stars correspond to the next 
generation of stars which have been previously enriched 
with metals. 

The computer program used for the calculations has 
been described by Klapp (1981, 1983) and Bahena 
(2006, 2007), with updated input physics. 

For the nuclear reaction rates we use the Nuclear 
Astrophysics Compilation of Reaction Rates (NACRE) 
by Angulo et al. (1999). A diffusion treatment for 
convection and semiconvection is used. For the opac- 
ity we have adopted the OPAL radiative opacitie s 



( Rogers and Iglesias 1119921 ; llglesias and Rogers 1 11993) 



3 Results 

3.1 Zero Age Main Sequence 

For our Pop I, II and III models, in Fig. [1] we show 
their ZAMS Hertzsprung-Russell (HR) diagram for the 
mass range 1 — 10, 000 A/©. The objective is to un- 
derstand the ZAMS structure differences as function of 
mass and metallicity. A large number of very detailed 
ZAMS models with different masses and metallicities 
have been calculated. 



The luminosity and effective temperature increases 
with mass. Low mass stars are located at the right lower 
part of the diagram, while massive and very massive 
stars are found in the left upper part, because they 
are the most luminous and hotter. Pop III stars are 
hotter than their Pop I and II counterparts, and so 
their locus on the HR diagram is shifted to the left 
upper part. Pregalactic stars are bluer than galactic 
stars. This is shown in Fig. 1 and also in Fig. 2 that is 
an amplification for the 100 — 10000Mq range. In Fig. 
3 we show a log g c - log T c diagram for a large mass 
and metallicity range. 

In Figs. 4 to 9 we show the main physical variables 
for the ZAMS models. All quantities are plotted as 
function of the mass. As the mass increases, ZAMS 
stars become bigger, brighter and less dense. With 
decreasing metallicity, Pop III stars get very hot and 
compact. All massive and very massive stars are dom- 
inated by radiation pressure and develop a large con- 
vective core. In all cases, however, their luminosity is 
below the Eddington upper luminosity limit. 

Then, we present ZAMS models in the metallicity 
range from Z = 10~ 10 to 10~ 6 for M = 100, 250, 500, 
750 and lOOOAf©, which are shown in Figs. 10 to 15. 
The central density increases with decreasing metallic- 
ity. The most massive stars are less dense. MS and 
VMS develop a large convective core, but its size de- 
creases slowly with decreasing metallicity. Stellar ra- 
dius is lower for low-metallicity. The radius also de- 
crease with decreasing mass, however, lower-mctallicity 
stars are more compact. Central temperature, and ef- 
fective temperature, increase with decreasing metallic- 
ity, and the most massive stars are hotter. Luminosity 
does not depend on metallicity but on the mass. The 
most massive stars are the most luminous stars. The 
Eddington luminosity factor T = L/L e dd, which is the 
ratio of the luminosity L to the Eddington luminosity 
L e dd does not depend on metallicity. The most massive 
stars are the ones closest to the Eddington luminosity, 
and the upper luminosity limit. 

3.2 Stellar structure 

We have calculated the evolution of 100, 150, 200, 250 
and 300M© stars with metallicities Z = 10~ 2 , 10 -3 , 
10 -6 , 10 -9 , and 10 -10 . Pop III stars have higher den- 
sity, temperature and pressure than their Pop I and II 
counterparts, and arc radiation pressure dominated and 
very luminous. 

Pregalactic stars have central temperatures of about 
log T c > 8 and effective temperatures of log T e jf < 5. 
A direct consequence of higher central temperatures is 
that they have higher energy generation rates. On the 
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Fig. 2 HR-diagram for ZAMS Pop I, II and III stars for the 
mass range 100 — 10, OOOA/© . Depending on metallicity the 
Pop III ZAMS is systematically shifted to higher effective 
temperature. 
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Fig. 3 ZAMS Pop I, II, and III Log £> c -Log T c plane for 
the mass range 1 — 10, 000M@. With increasing density and 
temperature, radiation pressure becomes more important. 
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Fig. 5 Ibidem. Effective temperature. 
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Fig. 6 Ibidem. Luminosity. 



Fig. 9 Ibidem. Gamma factor which is the ratio of the 
luminosity to the Eddington Luminosity. 
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Fig. 15 Ibidem. Gamma factor. 
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main sequence, lower metallicity stars produce slightly 
less energy. However, these stars are hotter than the 
others and so require higher temperatures to produce 
the same amount of energy. 

The most striking feature of the low metallicity stel- 
lar models is their atmospheric high temperature they 
are able to maintain. During hydrogen burning, these 
stars derive their nuclear energy from the inefficient pp- 
chains and the CNO-cycles. This is possible because a 
small fraction of carb on is produced during the pre- 



main sequence phase (jCastellani et al. 1 119831 ). Some 



carbon is also generated by the triple -a process be 



fore t he star reaches the main sequence (jBromm et al. 
2001 ). That is, in the absence of metals nuclear burn- 



ing proceeds in a non-standard way. First, the hydrogen 
burning occurs via the pp-chain. However, metal-free 
stars are hotter and very luminous reaching high cen- 
tral temperatures which are high enough for the simul- 
taneous occurrence of helium burning via the triple-a 
reaction. After a brief initial period of this process, 
a trace amount of heavy metals are formed and this 
makes that most of the energy generation rate during 
hydrogen burning comes from the CNO-cycles. 

Lower metallicity stars have a higher central tem- 
perature to support the star against gravity. This is a 
consequence of their high central density and tempera- 
ture, and because this type of stars are very compact. 



These tables summarize the most representative data 
of the models, and they include hydrogen and helium 
burning from the initial chemical composition, in 10% 
intervals of the helium mass fraction, to the end of he- 
lium burning when the helium mass fraction is approx- 
imately 0.01. 

VMS are hot and luminous and so are located in 
the left upper part of the Hcrtzprung-Russcll (HR) di- 
agram. Because Pop III stars are hotter than their 
enriched counterparts, their locus in the HR-diagram 
is shifted to the left upper part; pregalactic stars are 
bluer than galactic ones. 

For a given stellar mass, the evolution of a mas- 
sive star, i.e., its location in the HR-diagram, depends 
strongly on metallicity. Metal-free stars have unique 
physical characteristics and they exhibit high effective 
temperatures and small radii. In relationship with their 
cosmological consequences, metal-free models are im- 
portant for predicting the ionizing photon production 
of the first generation of stars. 

Pregalactic stars were denser and hotter than galac- 
tic Pop III stars. Some of their main physical vari- 
ables during the hydrogen and helium burning phases 
are shown in Figs [15] to 21. 

The following are the main properties of the studied 
evolutionary models for galactic and pregalactic very 
massive Pop III stars. 



3.3 Stellar evolution without mass loss 
3.3.1 Physical variables 

For the present work stellar evolution models without 
mass loss have been computed during the hydrogen and 
helium burning phases, for very massive galactic and 
pregalactic Pop III 100, 150, 200, 250 and 300M© stars 
with metallicities Z = 10~ 6 and 10 -9 , respectively. 

As an example of our evolution models, we list in 
Tables 1 and 2 some properties of 200M© galactic and 
pregalactic Pop III stars. 

The first two columns give the lifetime tq, given in 
units of 10 6 years, and the helium mass fraction Xn e 
during the hydrogen and helium burning phases. Sev- 
eral physical variables in logarithmic units are included 
such as the central density p c , the central temperature 
T Cl the luminosity L (in solar units), the effective tem- 
perature T g, the radius R, and the nuclear energy gen- 
eration rate e nuc . Then, the other quantities listed are 
the convective core size q cc , which is the ratio of the 
mass of the convective core to the total mass of the 
star, the radiation factor (3 C in the centre of the star, 
defined by 1 — f3 = P r /P, where P r is the radiation 
pressure and P the total pressure, and the Eddington 
luminosity factor T. 



a) Central density 

Pop III stars with Z = 10~ 6 settle down on the main 
sequence with central densities of p c — 7.46, 5.97 and 
3.94 g cm" 3 for 100, 200 and 300M Q , respectively. For 
the same stellar masses, and Z = 10~ 9 , central densi- 
ties are p c = 32.18, 29.03 and 24.83 g cm" 3 , respec- 
tively. That is, the central density is higher for lower 
metallicity stars, and so, Pop III pregalactic stars are 
denser than galactic stars. However, the central density 
decreases with increasing stellar mass. 

During hydrogen burning, the central density in- 
creases slowly but then by the end of this burning phase 
its increase rate grow significantly. The central density 
increases during the transition from hydrogen to helium 
burning and keeps increasing during the whole helium 
burning phase. 

b) Central temperature 

The VMS high energy requirements, demand high cen- 
tral temperatures in order to be able to maintain their 
structure and energy output. That is, the central 
temperature is high for more massive stars, but with 
lower metallicity the stars have an even higher tem- 
perature. At the beginning of hydrogen burning, the 



8 



Table 1 Physical variables and quantities for 200Mq galactic Pop III stars with initial metallicity Z — 10 , without 
mass loss and with no-rotation, during the hydrogen and helium burning phases. 
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Table 2 Physical variables and quantities for 200 Mq pregalactic Pop III stars with initial metallicity Z = 10 , without 
mass loss and with no-rotation, during the hydrogen and helium burning phases. 
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Fig. 16 Central density for 100M© (solid line), 150M© 
(dash-dot-dot-dot-dash), 200M© (dashes), 250M Q (dash- 
dot-dash) and 300M© (dots) pregalactic Pop III stars with 
metallicity Z = 10~ 9 , without mass loss, during the hydro- 
gen burning. 



Fig. 19 Central density for 100M Q (solid line), 150M Q 
(dash-dot-dot-dot-dash), 200M© (dashes), 250M© (dash- 
dot-dash) and 3OOM0 (dots) pregalactic Pop III stars with 
metallicity Z — 1CP 9 , without mass loss, during the helium 
burning phase. 
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Fig. 17 Ibidem. Effective temperature. 



Fig. 20 Ibidem. Effective temperature. 
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Fig. 18 Ibidem. Radius. 



Fig. 21 Ibidem. Radius. 
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central temperature with Z = 1CT 6 is T 7 = 7.02, 7.51 
and 7.76, for 100, 200 and 3OOM stars, respectively, 
where T n = T/10™ K. For Z = 10~ 9 , the corresponding 
temperatures are T 8 = 1.14, 1.30 and 1.36, respectively. 
During hydrogen burning, T c slowly increases until the 
end of this burning phase when it begins to increase 
to high values and continues to do so during the whole 
helium burning phase. 

c) Nuclear energy generation 

During hydrogen burning the main energy source is 
given by the CNO-cycles because of their strong tem- 
perature dependence. In lower metallicity stars these 
cycles are activated just after a short helium burn- 
ing phase that produces enough CNO elements for the 
CNO-cycles to operate. 

With increasing stellar mass, the central tempera- 
ture is higher and so the nuclear energy generation. 
However, for pregalactic stars it increases with decreas- 
ing metallicity. As a consequence of the larger cen- 
tral temperatures, CNO-cycles carbon production is en- 
hanced at an earlier evolutionary phase than in less 
massive stars. The nuclear energy generation increases 
during hydrogen burning. But, for galactic Pop III stars 
with metallicity Z — 10 -6 , the energy generation rate 
decreases when hydrogen tends to be exhausted. When 
the existing hydrogen mass fraction is about ~ 0.02, 
and the helium mass fraction reaches ~ 0.98, the nu- 
clear energy generation decreases. For pregalactic stars 
with Z = 10~ 9 , during the transition from hydrogen 
to helium burning there is only a small decrease in the 
nuclear energy generation rate, i.e. the transition from 
hydrogen to helium burning is smooth, quite different 
than for galactic stars. 

There arc different nuclear generation rates at the 
end of hydrogen burning because the central temper- 
ature in galactic Pop III stars is not high enough for 
helium ignition. This situation does not take place 
in the case of pregalactic stars precisely because they 
have a higher central temperature at the end of hydro- 
gen burning. Then, the transition to the next burning 
phase occurs in a smoother form. In the other case, 
a strong explosive helium flash takes place during the 
transition to helium burning. At this moment, the star 
is contracting, consequently, it is heated and reaches an 
appropriate central temperature to ignite helium. 

When helium is ignited in the core of the star, the 
nuclear energy generation rate increases rapidly dur- 
ing the transition to helium burning, then reaches a 
maximum and decreases towards the end of this burn- 
ing phase. This occurs for galactic stars with metal- 
licity Z = 10~ 6 . In the case of pregalactic stars, with 



metallicity Z = 10~ 9 , this transition takes place very 
smoothly because the stars are hot enough to ignite 
helium immediately after exhausting hydrogen. Then 
the generation rate increases towards a maximum and 
decreases by the end of helium burning. 

d) Luminosity 

Very massive Pop III stars are very luminous. For 100, 
200 and 300 Af© with metallicity Z = 10~ 6 , at the 
beginning of the main sequence, their luminosity are 
L = 1.22xl0 6 , 3.59xl0 6 and 6.31 xlO 6 L , respectively. 
For stars with Z = 10~ 9 , the corresponding figures are 
L = 1.27 x 10 6 , 3.71 x 10 6 and 6.51 x 10 6 L G . Dur- 
ing hydrogen burning the stellar luminosity increases 
slightly and varies smoothly during the transition from 
hydrogen to helium burning. Then it remains practi- 
cally constant during helium burning. 

e) Effective temperature 

Pregalactic Pop III stars are hotter than galactic stars. 
At the beginning of the main sequence, galactic stars 
with Z — 10~ 6 have an effective temperature T g = 
69582, 75586 and 77988 K for 100, 200 and 3OOM , 
respectively. For the same stellar masses, pregalactic 
stars with Z = 10~ 9 have T eS = 89170, 100090 and 
103390 K, respectively. These different values are due 
to their higher central temperatures and the different 
mechanisms to drive nuclear burning. 

The effective temperature continuously decrease dur- 
ing hydrogen burning, until the transition to helium 
burning that it increases. Then, the effective tempera- 
ture start to decrease but then remains practically con- 
stant until the end of this burning phase when the ef- 
fective temperature decreases even more. 

f) Radius 

Very massive Pop III stars are compact. On the main 
sequence, galactic stars with metallicity Z = 10 -6 have 
initial radii R = 7.67, 11.15 and 13.88i? for 100, 200 
and 300Mq, respectively. For pregalactic stars with 
Z = 10~ 9 , their radii are R = 4.75, 6.46 and 8.02i? o , 
respectively. 

As the effective temperature decreases, the radii in- 
creases. At the end of hydrogen burning, for galac- 
tic 100, 200 and 3OOM stars with initial metallicity 
Z = 10~ 6 their radii are R = 20.82, 33.40 and 43.39i? , 
respectively; and for pregalactic stars with initial metal- 
licity Z = 10~ 9 , R = 12.62, 18.61 and 24.18i? , respec- 
tively. 

On the contrary to effective temperatures, during 
the transition from hydrogen to helium burning, radii 
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slowly increase and remain almost constant during he- 
lium burning until the end when they increase again. 
This is, prcgalactic lower metallicity stars are more 
compact because they are hotter than their galactic 
counterparts. 

g) Convective core 

A convective core is always present from the ZAMS to 
the end of the helium burning phase. Very massive Pop 
III stars develop a large convective core. For 100, 200 
and 300Mq stars with Z = 10~ 6 , the convective core 
size at the beginning of the main sequence is q cc = 0.81, 
0.89 and 0.92, respectively. For the same stellar masses 
and Z = 10~ 9 , q cc = 0.79, 0.88 and 0.91, respectively. 
The core is larger for higher metallicity. 

The convective core size is larger for higher masses 
and increases during hydrogen burning. In fact, the 
studied stars are almost fully convective during hydro- 
gen burning. At the end of this burning phase the 
stars contract while forming a helium core. For the 
masses above mentioned and metallicity Z = 10~ 6 , at 
the end of hydrogen burning, the convective core size 
q cc = 0.42, 0.47 and 0.48 for 100, 200 and 3OOM , re- 
spectively. For Z = 10~ 9 , q cc = 0.40, 0.42 and 0.44, 
respectively. Then, stars with Z = 10~ 6 form a helium 
core of A/hc = 42.5, 98.3 and 142.7M©; with Z = 10" 9 , 
the core masses are Mu e = 39.7, 93.5 and 142. 4M Q , 
respectively. 

Because for galactic stars with metallicity Z = 10~ 6 
the transition from hydrogen to helium burning is ex- 
plosive, they suddenly contract, affecting momentarily 
their core size which rapidly increases, but then de- 
creases while forming a helium core. That is, the ex- 
haustion of hydrogen in the centre of the star causes 
a progressive contraction of the star and the shrink- 
ing of the convective core which finally vanishes when 
X c <~ 10~ 3 . During the transition, the energy released 
by the star is supplied by the gravitational contrac- 
tion. In the pregalactic case with Z = 10 -9 the burning 
transition is very smooth and the stars contract imme- 
diately forming a helium core. 

During helium burning, for 100, 200 and 300M Q 
stars with Z = 10~ 6 , the convective core is q cc ~ 0.40, 
0.45 and 0.47, respectively. For Z = 10~ 9 , g cc - 0.39, 
0.43 and 0.45, respectively. That is, a carbon core 
mass of M c = 40.5, 90.9 and 13O.6M , respectively, 
is formed for stars with Z = 10~ 6 , while for Z = 10 -9 , 
M c = 39.1, 86.4 and 134.7M Q , respectively. 

h) Radiation pressure 

Very massive Pop III stars are dominated by radi- 
ation pressure. At the centre of the stars we have 



1 - p c = 0.45, 0.58 and 0.64 for 100, 200 and 300Af Q 
galactic stars with metallicity Z = 10 -6 . That is, the 
contribution of radiation pressure to the total pressure 
is considerable. For pregalactic stars with metallicity 
Z = 10~ 9 we have that 1 - f3 c = 0.45, 0.58 and 0.64, 
respectively. 

The contribution of the central radiation pressure in- 
creases during hydrogen burning. After the transition 
from hydrogen to helium burning, the central radia- 
tion pressure contribution increases and then remains 
almost constant. 

i) The Eddington luminosity 

Very massive Pop III stars evolve during hydrogen 
burning below the Eddington upper luminosity limit. 
For Z = 10 -6 , at the beginning of hydrogen burn- 
ing, the ratio T = 0.60, 0.73 and 0.79 for 100, 200 and 
300M Q , respectively; and for Z = 10~ 9 T = 0.60, 0.74 
and 0.80, respectively. 

During hydrogen burning, this ratio increases and 
has a maximum close to the end of this burning phase, 
approximately when the helium mass fraction is 0.98 
and then decreases slightly as the star reaches helium 
ignition. After the transition from hydrogen to helium 
burning, T increases when the helium mass fraction is 
approximately 0.98. Furthermore, this ratio decreases 
slightly but then remains almost constant during the 
helium burning phase. 

3.3.2 Evolutionary tracks 

Fig. HU shows evolutionary tracks in the HR-diagram 
for 100, 150, 200, 250 and 300M o galactic and pregalac- 
tic Pop III stars with metallicity Z = 10 -6 and 10 -9 , 
respectively. The most massive stars are the hotter and 
most luminous. Luminosity decreases with lower metal- 
licity. For the same stellar mass, galactic and pregalac- 
tic Pop III stars have similar luminosities. All stars 
settle down on the main-sequence with a high effective 
temperature and luminosity. During hydrogen burn- 
ing all stars increase their luminosity while the helium 
mass fraction increase and both the central density and 
temperature increase. 

For stars with metallicity Z = 10~ 6 , the transi- 
tion between nuclear burnings is explosive while the 
star contracts. At this moment, the effective tempera- 
ture and luminosity increase. Then, the luminosity re- 
mains almost constant while the effective temperature 
decreases. However, for different stars their luminos- 
ity and effective temperature decrease with decreasing 
mass and metallicity. For pregalactic stars the hydro- 
gen to helium burning transition occurs very smoothly 
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Fig. 22 Evolutionary tracks in the HR-diagram for 100, 150, 200, 250 and 300Af© Pop III stars with metallicity Z = 10" 
and Z = 10 -9 , respectively and without mass loss during the hydrogen and helium burning phases. 



because they are hotter and their central temperature 
is high enough to ignite helium promptly. 

Galactic and prcgalactic stars evolve with differ- 
ent lifetimes. For 100, 200 and 3OOM galactic Pop 
III stars, their lifetimes during hydrogen burning are 
3.07566, 2.42932 and 2.20822 megayears, respectively. 
For prcgalactic stars, their lifetimes for the same masses 
are 2.87012, 2.28397 and 2.09287 megayears, respec- 
tively. During helium burning, stellar lifetimes arc 
shorter than during hydrogen burning. 

The initial hydrogen and helium burning phases take 
place at the blue side of the HR-diagram. When the 
stars evolve they move toward the red. However, the 
presently studied galactic and pregalactic very mas- 
sive Pop III stars, evolving without mass loss and no- 
rotation, do not experience the asymptotic giant branch 
(AGB) phase. The most massive zero-metal stars tend 
to be cooler but the temperature remains above ~ 10 4 
K. 

The HR-diagram evolution of massive lower metallic- 
ity Pop III stars is different than for Pop I and II stars, 
because they do not evolve to red-giants before core 



collapse (lUmeda et al. I|2000D. In the q uoted range of 



masses used by ( Castellani et al. 1.983) 



massive zero- 



metal stars fail to reach the red giant region either at 
the hydrogen or helium burning phases. 



In the Hcrtzprung- Russell diagram the locus of very 
massive Pop III stars is in the left upper part. These 
stars are hotter and very luminous. Pregalactic stars 
are hotter than galactic Pop III stars. Then, stars with 
lower metallicity are shifted to the left because they are 
bluer than the others. 

k) p — T plane 

The p — T plane describes the state of the gas in the in- 
nermost stellar regions and it is important for the diag- 
nosis of the stellar structure and evolution. The evolu- 
tion of the central conditions determine the boundaries 
in which the equation of state is dominated by different 
pressure components, i.e., ideal gas, radiation pressure 
or electron pressure. 

According to the zones of the equation of state of 
an electron gas, the studied stars occupy the upper 
loci of a non-degenerate and non-rclativistic gas. In 
this zone there is a boundary at which pair production 
could become important. Regarding the gas in ther- 
modynamic equilibrium, these stars are dominated by 
radiation pressure. 

The central conditions varies in the range log p c <~ 
1 — 3 and log T c = 8.0 — 8.6 for galactic stars with 
Z = 10~ 6 . For pregalactic stars with Z = 10 -9 , both 
the central temperature and density are higher. The 
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evolution of the central conditions can be described by 

1/3 . 
T c ~ p c ■ This behaviour depends on the equation of 

state and it is similar for models of any metallicity. 

We distinguish two regions in the p — T plane in 
correspondence to the stages of gravitational contrac- 
tion of stellar cores between nuclear burnings. In the 
second region, corresponding to helium burning, the 
central density and temperature increases more than 
during hydrogen burning. 



4 Discussion 

The evolutionary tracks of our massive Pop III stars are 
shifted to the left upper part of the HR-diagram as in 
Tumlinson et al. (2003) models. Then, stars evolve to 
the red with increasing luminosity and decreasing effec- 
tive temperature. In both cases, luminosities are simi- 
lar but, in the present case, final effective temperatures 
are slightly higher. This is probably due to the differ- 
ent parameters used, chemical composition and specific 
implementation of physical processes, e.g., convection. 

a) Nuclear lifetimes 

Our models are hot and luminous and have short life- 
times. For 100 M Q stars with Z = 10~ 6 and 10~ 9 , their 
lifetimes during hydrogen burning are th = 3.09368 
and 2.92147 megayears, respectively. Lifetimes during 
helium burning are of the order of ~ 10% of their life- 
time during hydrogen burning. Our results for galactic 
and prcgalactic hydrogen burning Pop III stars are in 
good agreement with other authors. 

For M > 300M Q , even metal- free stars evolve toward 
T c ff < 10 K and eventually become red supergiants, as 
the hydrogen burning s hell becomes more ac tive with 
increasing stellar mass (jBaraffe et al. 1 120011 ) . In the 



present work, the most massive metal-free stars tend 
to be cooler and likely could become red supergiants. 
However, for a helium mass fraction equal to ~ 0.01, 
they maintain high effective temperatures. In all cases, 
stars do not reach low effective temperatures. Evolving 
stars without mass loss and with no-rotation, fail to 
reach the AGB phase. At the end of helium burning, 
the studied galactic stars are hotter than log T e g ~ 4.0 
and pregalactic stars are hotter than T e s ~ 4.6. 

b) Nuclear energy generation 

The peculiar behaviour of low metallicity stars was first 
pointed out by Ezer (1961) and their structure during 
hydrogen and helium burning has been investigated by 
several authors. 



According to the results presented here the stars be- 
gin to settle down on the main-sequence with higher 
initial central temperatures of order log T c ~ 8. Then, 
the onset of the 3a reaction occurs at earlier stages of 
hydrogen burning. This is because the 3a reaction re- 
quires a much higher temperature than the pp-chains 
and ignites at the beginning of the hydrogen burning 
phase. In the present models, after a brief initial period 
of 3a burning, a trace amount of heavy elements has 
been formed. Then, the stars expand and follow the 
CNO-cycles. 

In Fig. [23] the nuclear energy generation as function 
of the temperature is shown, and Fig. [24] shows this 
energy generation during hydrogen burning for both 
galactic and pregalactic Pop III stars. 

c) Convective core size 

In their models, Marigo et al. (2001) found that if 
the convective core grows, it eventually reaches the H- 
shcll and engulfs some hydrogen-rich material, which is 
rapidly burnt via the CNO-cycles. This causes a flash 
that expands the core, so that central helium burning 
weakens and the convective core recedes temporarily 
(in mass). After the flash has occurred, the convective 
core starts growing again. 

In the present models, the same picture takes place 
for galactic Pop III stars but it does not in the pre- 
galactic case. In this case, the central temperature at 
the transition of nuclear burnings is higher than in the 
first one. 

Hydrogen exhaustion in the centre of the star causes 
a progressive contraction of the star and the shrink- 
ing of the convective core which finally vanishes when 
X c ~ 10~ 3 , so that no standard over all contraction 
phase is found ([Castellani et al. 1 1 1983! ) . A maximum 



in the energy released by gravitation occurs during this 
burning phase, when X c = 10 -9 , and 59% of the energy 
released by the star is supplied by contraction. Ignition 
of the full triple-a cycle is only slightly delayed with re- 
spect to hydrogen burning ignition in a shell, and once 
again the gravitational contraction supplies energy to 
the star. Once the 3a chain has become fully efficient, 
a convective core is again developed. 

For the present work, pregalactic Pop III stars have 
higher temperatures than galactic stars allowing a soon 
ignition of the 3a reactions. Then the transition from 
hydrogen to helium burning is smooth because the con- 
vective core does not vanish. 

d) Eddington luminosity 

Stellar models may become gravitationally unbound, 
i.e., the rate of energy outflow at the surface exceeds 
the corresponding Eddington luminosity. 
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Fig. 24 Main sequence nuclear energy generation rate as function of the hydrogen mass fraction for 2OOM0 galactic (left 
panel), and pregalactic (right panel) Pop III stars, with metallicity Z = 10 -6 and Z = 10 -9 , respectively. 
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In the present models, the Eddington T factor is 
calculated by estimating the ratio between the cur- 
rent stellar luminosity and the Eddington luminosity. 
In both the hydrogen and helium burning phases, the 
studied cases of 100, 200 and 300Af© galactic and pre- 
galactic Pop III stars evolve below the upper Eddington 
luminosity limit. During hydrogen burning, the T fac- 
tor increases, reaches a maximum toward the end of this 
burning phase and then decreases. The most massive 
stars evolve during the helium burning phase closest to 
the Eddington limit. 



5 Conclusions 

For the present work a large number of ZAMS models 
have been calculated showing the main physical vari- 
ables as function of the mass, and the metallicity. As 
the mass increases, ZAMS stars get bigger, brighter, 
and less dense. Pop III stars get very hot and compact 
as metallicity decreases. 

When metal-free stars settles down on the main- 
sequence, they have smaller radii, hotter cores, and 
higher effective temperatures than metal-enriched stars. 
Like their counterparts, lower metallicity stars become 
systematically cooler, larger and more luminous over 
their hydrogen burning lifetimes. The most massive 
stars have shorter lifetimes than less massive stars. 

To study their properties and internal structure, stel- 
lar structure models on the main sequence have been 
calculated, emphasizing the case of metal-free stars 
which are compared with their Pop I and II coun- 
terparts. Pop III stars are more centrally condensed, 
denser and hotter. On the main sequence their different 
internal regions are always below the upper Eddington 
luminosity limit. 

The most important feature of the metal-free models 
is the high temperature they maintain in their photo- 
sphere. These stars have a high ionizing photon pro- 
duction rate. The ionization caused by these stars is a 
direct result of their high effective temperatures. For 
M > 100M Q , the studied metal- free stars have effec- 
tive temperatures T e g ~ 10 5 K. Consequently, they are 
very efficient at producing photons capable of ionizing 
hydrogen and helium. According to Bromm, Kudritzki 
and Loeb (2001), very massive Pop III stars might have 
played a significant role in the IGM rcionization of hy- 
drogen and helium at high redshifts. 

Nuclear burning proceeds in a non-standard way. 
For stars with no metals, the pp-chains and the CNO- 
cycles do not provide enough energy to support the star, 
so when it reaches the main sequence keeps contract- 
ing until the 3a process starts to operate, halts the 



contraction and in a very small time produces enough 
CNO elements so that the energy from the CNO-cycle 
is capable of supporting the star, which expands and 
settles down to the main sequence. The stars maintain 
core temperatures in excess of T c > 10 8 K which are 
high enough for the simultaneous occurrence of the pp- 
chains, the CNO-cycles and helium burning via the 3a 
process. 

For these stars, the radiative opacity in their en- 
velopes is reduced, and their core temperature is high, 
then the first stars are hotter and smaller than metal- 
enriched stars. On the other hand, the opacity of stellar 
matter is reduced at low metallicity, permitting steeper 
temperature gradients and more compact configura- 
tions at the same mass. 

Very massive (M > 100M Q ) galactic and prc- 
galactic Pop III stars develop large convective cores 
with important helium core masses Mh ~ 40 Af©. 
These quantities are impo rtant for explosion calcula- 
tions (jUmeda et al. 1120001) and synth etic derivation of 



the SN yields (jPortinari et al. 1119981) . Semiconvcction 



does not greatly affect the stellar structure during the 
main-sequence phase. However, during the shell hy- 
drogen burning and helium burning phases, it plays a 
significant role, and the evolutionary results depends 
on the adopted criterion and the input physics of the 



models ( Chiosi and Maeder 19861) 



We have calculated evolutionary models for lower 
metallicity very massive Pop III stars. The evolution 
of these stars is similar to that of metal enriched stars 
but now the evolutionary tracks are shifted to the left 
of the HR-diagram, i.e. the models are bluer than their 
metal enriched counterparts. This is because massive 
Pop III stars are hotter and very luminous. Pregalacic 
stars are denser and hotter than galactic Pop III stars. 
However, during their evolution these stars are more 
luminous than the first ones. 

During the hydrogen burning phase, very massive 
galactic and prcgalactic Pop III stars evolve almost fully 
convective. They form a large core with a different 
structure depending on metallicity. Very massive stars 
are dominated by radiation pressure and electron scat- 
tering. As a consequence of the high radiation pressure, 
the convective core tends to be larger at higher masses 
and to expand as the star evolves. 

Because the stars radiate near the Eddington limit, 
radiation pressure due to electron scattering opacity 
can become substantial. For pregalactic Pop III stars, 
with the metallicity Z = 10~ 9 considered in this work, 
they develop a helium core Mn e = 39.7, 93.5 and 
142. 4Mq, corresponding to initial stellar masses of 100, 
200 and 300 A/©, respectively. Then, according to Fryer 
et al. (2001) and Heger and Woosley (2002), these stars 
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will likely explode by pair-instability supernovae for the 
200 and 300M© cases or collapse into black holes for 
the IOOM0 case. However, M < 13OM0 stars could 
explode like hypcrnovac. 

The present evolutionary models have been calcu- 
lated without mass loss and with ro-rotation, in accor- 
dance with other existing points of view. However, the 
uncertain role of mass loss can be viewed as one of the 
major systematic uncertainties remaining in the study 
of metal-free stars. In subsequent papers we will extend 
the present study to mass losing and rotating models. 
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